SCIENTIFIC 

REPORTS 




SUBJECT AREAS: 

BACTERIA 

ENVIRONMENTAL 
MICROBIOLOGY 

GENE EXPRESSION 

PROTEOMICS 



Received 
29 March 201 1 

Accepted 
16 June 2011 

Published 
29 June 2011 



Correspondence and 
requests for materials 
should be addressed to 
JJ.P. (jacob.parnell® 
usu.edu] 



Time-course analysis of the Shewonello 
omozonensis SB2B proteome in response 
to sodium chloride shock 

J. Jacob Parnell 1 , Stephen J. Callister 2 , Giovanni Rompato 1 , Carrie D. Nicora 2 , Ljiljana Pasa-Tolic 5 , 
Ashley Williamson 3 & Michael E. Pfrender 4 

'Center for Integrated BioSystems and Department of Biology, Utah State University, Logan, UT 84322, 2 Biological Sciences 
Division and the Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA 99352, 
3 Logan High School, Logan, UT 8432 1 , 4 Department of Biological Sciences, University of Notre Dame, Notre Dame, IN 46556. 

Shewanellae are microbial models for environmental stress response; however, the sequential expression of 
mechanisms in response to stress is poorly understood. Here we experimentally determine the response 
mechanisms of Shewanella amazonensis SB2B during sodium chloride stress using a novel liquid 
chromatography and accurate mass-time tag mass spectrometry time-course proteomics approach. The 
response of SB2B involves an orchestrated sequence of events comprising increased signal transduction 
associated with motility and restricted growth. Following a metabolic shift to branched chain amino acid 
degradation, motility and cellular replication proteins return to pre-perturbed levels. Although sodium 
chloride stress is associated with a change in the membrane fatty acid composition in other organisms, this is 
not the case for SB2B as fatty acid degradation pathways are not expressed and no change in the fatty acid 
profile is observed. These findings suggest that shifts in membrane composition may be an indirect 
physiological response to high NaCl stress. 

Shewanellae are models for studying the genomic and phenotypic response to environmental stress because 
of their ability to inhabit and adapt to a wide range of environmental conditions 1 . Shewanellae are capable of 
both aerobic and anaerobic respiration utilizing diverse electron acceptors (fumarate, thiosulfate, nitrite, 
nitrate, iron, chromium, manganese, and uranium) 2 6 . Consequently, Shewanellae have robust sensing and 
regulatory systems that allow rapid response and acclimation to changing environments 2 . Their versatility and 
plasticity provides an opportunity to investigate how microorganisms respond to changes in typical envir- 
onmental conditions, such as temperature, dessication, and osmotic stress 7 9 . Although they are typically mar- 
ine-associated, Shewanella have been detected in environments ranging from fresh water 6, 10 to hypersaline 
environments 1113 and tolerate a wide range of salt concentrations 9 . 

The genotypic and phenotypic response of Shewanella amazonensis SB2B (SB2B) to osmotic stress has not 
been described. This bacterium was isolated from shallow-water marine deposits derived largely from the 
Amazon River delta 14 . The physical mixing of these deposits by wave action combined with pore water transport 
exposes these deposits to periodic fluctuations in natural salinity 14 . Environmental fluctuations in salinity of its 
native habitat suggest that S. amazonensis is able to mitigate the effects of salt toxicity and hyperosmolarity. To 
experimentally determine the response mechanisms of S. amazonensis SB2B, we analyzed the dynamics of the 
proteome from initial response to acclimation (resumption of growth) during sodium chloride stress. The 
genome of SB2B is sequenced, providing peptide sequence information for conducting high throughput pro- 
teomics analyses using the accurate mass and time (AMT) tags approach 1517 . We use novel approaches in time- 
course proteomics to gain insight to the sequence of events and coordination of mechanisms involved in salt 
tolerance such as initial signal transduction and regulation, turgor maintenance, motility, and metabolism. 

Bacterial genomic, and proteomic responses to sudden changes in osmolarity play an important role in the 
biology and ecology of microorganisms because hyper- and hypo- osmolarity are stresses most commonly 
encountered in natural environments 1819 . Microbes colonizing almost any environment face fluctuations in 
osmolarity 19 , and must be equipped with response mechanisms in order to survive the initial change as well as 
prolonged exposure to osmotic stress. The fluctuations in salinity are a key determinant in the composition of 
microbial communities in estuarine and brackish water ecosystems 20 such as those from whence SB2B was 
isolated 14 . Bacterial response to osmotic stress is a complex, distinct sequence of cellular events 21 and elucidation 



SCIENTIFIC REPORTS | 1 : 25 | DOI: 1 0. 1 038/srep00025 



1 



of dynamic cellular responses to salt stress benefit from the use of 
global techniques that allow high throughput analysis of proteins 
involved in the process. Proteomic methods are ideal for the analysis 
of global changes and have widely been used to enhance the know- 
ledge of differentially expressed functional and regulatory proteins. 
However, until recently, proteomic technology and analytical proto- 
cols lacked the quantitative foundation to provide a platform for 
performing detailed time-course studies 22 23 . The aim of this study 
is to characterize the acclimation of Shewanella amazonensis str. 
SB2B to high sodium chloride stress using novel time-course proteo- 
mic techniques. 



Sama.U-ll 



Results 

To determine SB2B's salt tolerance, we exposed cultures to concen- 
trations of salt ranging from 0.3M to 1.37M. The growth rate of 
SB2B was inhibited by approximately 50% at a salt concentration 
of 0.85M in growth cultures (Fig. la). Upon exposure of SB2B to 
0.85M salt, growth is arrested between approximately 60 to 90 min 
following sodium chloride pulse (Fig. lb). Accordingly, the 
response frame used throughout this study ranged from 0 to 
90 min, and was broken down to measure the initial response (up 
to 15 min), the intermediate response (15 to 60 min), and the 
sequence of events in longer-term acclimation (60 to 90 min) of 
SB2B to salt stress. 

A total of 845 (out of 3,645 protein coding sequences; 23.2%) 
proteins were identified from SB2B as expressed throughout all five 
time points measured in this study. One hundred forty-three of the 
detected proteins demonstrated more than a two-fold change in 
expression (relative to time 0) for at least one subsequent time point. 
Principal component analysis of the expression patterns of all differ- 
entially expressed proteins suggest two separate responses to high 
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Figure 1 | The effect of NaCI on S. amazonensis SB2B. (A) Growth rates 
in doubling IT 1 (y- axis) with increased salt concentration (x- axis). (B) 
Time-course plate count results of following NaCI pulse to 0.85M 
concentration. 



Principal Component 1 (60°o) 

Figure 2 | Principal component analysis of all differentially expressed 
proteins. Hierarchical clustered groups involved in branched chain amino 
acid degradation (green) and replication and division (red) account for 
most of the variability in expression patterns. 

salt (Fig. 2). Most of the variation within expression is described by 
proteins corresponding to arrested growth of SB2B and a shift to 
branched chain amino acid degradation. 

Within the first 15 min following 0.85 M salt exposure, SB2B 
induces a number of genes involved in protein turnover and tran- 
scription. The first 15 min of arrested growth of SB2B correlates to a 
significant decrease in the abundance of DNA polymerase 
(Samal310), which was also confirmed by RT-qPCR (Table 2). A 
significant decrease in ribosomal proteins (Sama2548, 3033, and 
1999) detected, also suggests arrested cell growth and a reduction 
in translation, as illustrated under 'primary repressed' in figure 3 
(although the trend is similar for other ribosomal proteins, 
Sama0202-0239 show no significant change in expression). Other 
replication (TopA; Sama3041, MutS; Samal045, and DNA helicase; 
Sama3058) and cell division (FtsK; Samal772, FtsZ; Sama0358) pro- 
teins were also lower following salt addition. TopA and FtsK levels 
were also confirmed using RT-qPCR (Table 2). In addition, we found 
significant expression of a methyltransferase (thiopurine S-methyl- 
transferase - Sama0543). Reverse transcriptase-quantitative PCR 
confirms decreased expression of DNA replication and translation 
proteins that correspond to the diminished growth rate. One of the 
most highly differentially expressed proteins (8.4-fold increase) 
within the first 15 minutes of salt exposure is a diguanylate cyclase 
(Sama2418). 

In addition to a pause in cell division proteins, SB2B induces (or 
represses) a pulse of proteases (proteases; Sama2056 and Sama3074) 
and nucleases (endonucleases; Sama0514 and Sama0278, exonu- 
clease SbcC; Samal268, and ribonucleases; Samal933 and 
Sama3067: see Supplementary Table 1). Immediately following the 
NaCI pulse, expression of the protease/chaperone ClpA (Sama2056) 
is lower and remains so until growth resumes between 60 to 90 min- 
utes following the applied stress. 

Although some fatty acid degradation proteins (Sama2071; oxi- 
doreductase, and Sama2167; fatty acid oxidation alpha subunit) were 
expressed, we did not observe a change in the fatty acid content of 
SB2B following exposure to 0.85 M NaCI (described below). While 
proteins involved in branched-chain amino acid degradation 
increased in relative abundance 15 to 30 min following NaCI expo- 
sure, proteins involved in the synthesis of branched chain amino 
acids did not significantly change in relative abundance. Fatty acid 
analysis of SB2B throughout acclimation to high salt conditions 
showed no significant difference in branched chain fatty acid com- 
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Table 1 | Functional genes and corresponding polymerase chain reaction primers used for cDNA quantitation in this study. 



Forward primer 



Reverse primer 



3-ketoacyl-CoA thiolase Sama0031 

Hydroxymethylglutaryl-CoA lyase Samal 358 

3-methylcrotonoyl-CoA carboxylase Sama 1359 

Short chain enoyl-CoA hydratase Sama 1 378 

Acetyl-coenzyme A synthetase Sama2079 

2-methylcitrate synthase Sama3295 

Aconitase Sama3296 

Flagellar biosynthetic protein FlhB Sama2284 

Flagellar hook-length protein FliK Sama2292 

Flagellar hook protein FIgE Sama23 1 5 

Cell division protein FtsZ Sama0358 

DNA-directed DNA polymerase Sama 1310 

DNA topoisomerase IVsubunitA Sama3041 



Branched Chain Amino Acid Degradation 



5 '-ACAGTTAACCGCCTCTGTGG-3 ' 
5 '-TAAAGGCCCGTAAAATGGTG-3 ' 
5 '-TGCTTGACGTCGATATGAGC-3 ' 
5 'ACGTATCGAAGGCAATACCG-3 ' 

TCA cycle 
5 '-AACCCGGTGAAGTGTACTGG-3 ' 
5 '-AACCTGGAAACCGAACACAG-3 ' 
5 '-CTCGAGAAGGATGCCTTTG-3 ' 
Motility 

5 '-TTCTGTTGCTGTGCAGTTCC-3 ' 
5 '-TTGATGCAAACTCTGCTTGG-3 ' 
5 '-CTGGACTGGTTCTGTGCTGA-3 ' 
Replication and Division 
5 '-CGCCGAAGAGCCTGTATATC-3 ' 
5 '-AGCACTGTTCCGGTAACCAC-3 ' 
5 '-AACCGGTGCAAAGACGTATC-3 ' 



5 '-TCGACACCACCGATGATAAA-3 ' 
5 '-GCCGTAGGTGTCGTGAAAGT-3 ' 
5 '-CTGCCATTTCACCAAATCCT-3 ' 
5 '-TCAGCTCATTCAGGCTCTCA-3 ' 

5 '-TGGGGATACCCTCATGGATA-3 ' 
5'-ATCGTGGCTGAAACGATACC-3' 
5 '-AAGGCCTTTTTGGTCCAGTT-3 ' 

5 '-GAGCTGCTTCTTGTGGTTCC-3 ' 
5 '-CCTGATATCGGCTGGATGTT-3 ' 
5 '-CCCAGAGGTTCTGTGGTGAT-3 ' 

5 '-AGCTCGTTCTTTGGTGCAGT-3 ' 
5 '-TGTTTGGGTATTGGCTGACA-3 ' 
5 '-CGTACGGGCAGATTTCTTGT-3 ' 



position (Table 2). The only significant decline for fatty acids 
once SB2B was exposed to high salt concentrations was 16:0 and 
16:1. 

Three malonyl acyl carrier proteins were detected (Samal978, 
1982, and 2154) but none showed any difference in expression due 
to high salt. Sama 1362 encodes an isovaleryl-CoA dehydrogenase 
(E.C. 1.3.99.10) specific for 3-methylbutanoyl-CoA in the complete 
degradation of leucine. We also observed an increased abundance of 
two protein homologs (Samal 377 and Sama 1870) that are butyryl- 
CoA dehydrogenases (E.C. 1.3.99.2) specific for isobutyryl-CoA and 
2-mefhylbutyryl-CoA in the complete degradation of valine and 
isoleucine, respectively. The alternative to degradation of these 
amino acids at this junction in the pathway results in synthesis of 
branched chain fatty acids. From this point in the degradation of 
branched chain amino acids, each step is expressed in a similar 
fashion and to a similar extent (Fig. 3). Both Samal378 and 1379 
encode an enoyl-CoA hydratase (E.C. 4.2.1.17) that can be used in all 
three amino acid pathways. To complete the degradation of leucine 



to acetoacyl-CoA, we found expressed proteins for the 3-methylcro- 
tonyl-CoA carboxylase (Samal359; E.C. 6.4.1.4), hydroxymethylglu- 
taryl-CoA lyase (Samal358; E.C. 4.1.3.4), and two expressed 
homologs for 3-oxoacid CoA-transferase (Samal357 and Samal457; 
E.C. 2.8.3.5). For the continued degradation of valine, we found an 
expressed 3-hydroxyisobutyrate dehydrogenase (Samal380; E.C. 
1.1.1.31) and an NAD -dependent aldehyde dehydrogenase 
(Sama2648; E.C. 1.2.1.3) that can shunt to succinyl-CoA as well as 
a methylmalonate-semialdehyde dehydrogenase (Samal 376; E.C. 
1.2.1.27) that can shunt to propanoyl-CoA. For the continued degra- 
dation of isoleucine, an acetyl-CoA acyltransferase (Sama0031; E.C. 
2.3.1.16) was expressed that yields propanoyl-CoA. We found that 
aconitase (Sama3296; E.C. 4.2.1.3), 2-methylcitrate synthase 
(Sama3295; E.C. 2.3.3.5), and 2,3-dimethylmalate lyase (Sama3294; 
E.C. 4.1.3.30) are grouped together on the chromosome and were 
highly expressed between 15 and 30 minutes following salt addition. 
In addition, we found similar expression of an acetyl- Co A synthetase 
(Sama2079; E.C. 6.2.1.1). 



Table 2 | Validation of proteome data using quantitative PCR. 

RT- PCR Proteomics 
1 5 min 30 min 60 min 90 min 15 min 30 min 60 min 90 min 







Branched Chain Amino Acid Degradation 










3-ketoacyl-CoA thiolase 


Sama0031 


-0.4 


0.9 1.9 


3.3 


-0.5 


0.1 


1.4 


1.2 


Hydroxymethylglutaryl-CoA lyase 


Samal 358 


-1.7 


0.7 3.6 


1.8 


-2.1 


2.9 


3.9 


4.1 


3-methylcrotonoyl-CoA carboxylase 


Samal 359 


-0.5 


1.8 2.1 


1.8 


-0.3 


2.1 


2.2 


2.3 


Short chain enoyl-CoA hydratase 


Samal 378 




0.2 2.2 


1.3 


0.1 


0.5 


1 .3 


1.4 








TCA cycle 












Acetyl-coenzyme A synthetase 


Sama2079 


-0.1 


2.2 3.1 


2.8 


-0.4 


3.7 


4.8 


5.1 


2-methylcitrate synthase 


Sama3295 


-0.4 


1.1 -0.3 


1.6 


0.1 


1.8 


2.0 


1 .8 


Aconitase 


Sama3296 


1.7 


5.8 1.6 


3.2 


-0.2 


1.6 


1.8 


1.7 








Motility 












Flagellar biosynthetic protein FlhB 


Sama2284 


-0.3 


-0.5 -1.9 


-4.5 


nd 


nd 


nd 


nd 


Flagellar hook-length control protein 


FliK Sama2292 


-5.9 


-7.9 -1.0 


-3.9 


nd 


nd 


nd 


nd 


Flagellar hook protein FIgE 


Sama23 1 5 


0.1 


-3.0 -0.8 


-0.2 


nd 


nd 


nd 


nd 








Replication and Division 












Cell division protein FtsZ 


Sama0358 


-1.6 


-2.1 -0.5 


-0.6 


-0.2 


-0.2 


-0.4 


-0.2 


DNA-directed DNA polymerase 


Samal 3 1 0 


-2.3 


-5.3 -1.2 


-6.2 


-6.1 


1.0 


-7.0 


3.0 


DNA topoisomerase IV subunit A 


Sama3041 


-0.2 


-1.7 0.1 


-0.8 


-0.3 


-0.3 


-2.1 


-0.3 



*nd - Not determined 
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Figure 3 | Cellular model of differentially expressed pathways and proteins in response to NaCl stress in SB2B. Pathways and proteins are represented 
with expression values of 1 5 min, 30 min, 60 min, and 90 min following exposure to NaCl relative to no exposure as bars next to each respective protein. 
Immediate induction (primary induced) proteins on the left suggest a signal cascade that influences cellular motility. Immediate repression (primary 
repressed) proteins in the middle suggest the arrest of DNA replication and cell division. Both the primary induced and primary repressed proteins 
approach pre-stressed expression levels at 90 min after NaCl pulse. Secondary induced proteins (right) involved in branched chain amino acid 
degradation are induced between 1 5 and 30 min following addition of NaCl; their expression value remains high past the 90 min time frame of this study. 



By 90 min following salt exposure, resumption of growth in SB2B 
correlates with the abundance of proteins involved in DNA replica- 
tion and cell division to pre-osmotic shock levels. The methyltrans- 
ferase (Sama0543), diguanylate cyclase (Sama2418), endonuclease 
(Sama0278), and exonuclease SbcC (Samal268) that were expressed 
as described earlier, decreased in abundance well below pre-per- 
turbed levels. Conversely, DNA polymerase (Samal310) and the 
DNA mismatch repair protein MutS (Samal045) go from decreased 
abundance following exposure to salt (15 min) to >2-fold increase 
(90 min) relative to pre-salt exposure. 

Physical measurements of motility confirmed that SB2B decreases 
motility upon exposure to high salt concentrations (Fig. 4). However, 
the protein expression patterns for decreased motility appear late in 
the acclimation to salt stress. 

Discussion 

Time course proteome analyses of the response of SB2B suggest an 
immediate, but brief motile osmotactic response to salt stress. This 
motility response corresponds with a pause in growth and cell divi- 
sion observed in growth curves (Fig. 1) as well as in the low express- 
ion of replication and division proteins. The secondary response of 
SB2B is a shift to branched- chain amino acid metabolism prior to 
resuming growth. Figure 3 illustrates the proposed model of the 
primary response, and secondary response of SB2B to salt stress. 
Principal component analysis of all differentially expressed proteins 
shows that 90% of the variation in expression is accounted for by two 
components (60% for component 1, 30% for component 2; Fig. 2). 



Most of the proteins involved in the variation of expression fall into 
the two categories mentioned above (signal transduction and 
arrested replication/division, and branched chain amino acid degra- 
dation). 

Signal transduction systems play a major role in the acclimation of 
microorganisms to changes in osmolarity and histidine kinases have 
been identified as osmosensors in both prokaryotic 24 25 and eukar- 
yotic cells 26 . Sama2418 is a diguanylate cyclase response regulator 
containing a CheY-like receiver domain and a GGDEF domain that 
may affect the direction of rotation and of the flagellar motor to 
control the motile behavior of the bacterium 25 . In several species, 




1% 3% 
Agar Concentration 



Figure 4 | Physical measurements of SB2B motility on salt. Mean 
diameter (in mm; ± standard deviation) of SB2B growth on soft agar (1% 
agar in LB medium; left) and solid agar (3% agar in LB medium; right) at 
low salt concentration (0.2M; dark) and high salt concentration (1.0M; 
light). 
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Amidohydrolase Sama00S7 

DNA binding protein fis Sama3226 

TonBdcpendent receptor Sama0131 

ATP-dependentClp protease Sama2056 

Thiopurine S-methyltransfe rase Sama0543 

Lipid-A disaccharide synthase SamallS2 

16SrRNAmethyltransferase Sama254S 

Chorismatemutase SamaOS9S 

Uncharacterized protein Sama2617 

Diguanylate cydase Sama241S 

tRNA modification GTPase MnmE SamaOOOS 

Serine threonine kinase SamaI340 

Endonudease Sama027S 

Uncharacterized protein Samal5S2 

Small protein A Sama23$5 

DNA mismatch repair MutS Samal045 

Thymidylate kinase Sama2043 

Response regulator receiver protein Sama2S0O 
Ribosome small subunit-dependent GTPase A Sama3033 

DNA-directed DNA polymerase Samal 3 10 

FAD-binding oxidoreductase Sama2071 

RNAse R Sama3067 

Uncharacterized protein Sama2216 

AcetyiCoA synthetase Sama20~9 

HydroxymetbylglutarylCoA lyase Samal35S 

4-aminobuty rate aminotransferase Sama2636 

Prolyl oligopeptidase family Sama202S 

J-methylcrotonoyl-CoA carboxylase Samal359 

2- methylcitrate synthase Sama3295 
2.3-dimethvlmalate lyase Sama3294 
Methylmalate semialdehyde dehydrogenase Samal3"6 

3- oxoacid-CoA transferase Saraal35~ 
Aconitase Sama3296 
PropionyfCoA carboxylase Samal361 



Figure 5 | Hierarchical clustering of protein expression profiles. All differentially expressed proteins (143) clustered (left). Closely clustered groups 
containing replication/division proteins (top right) and branched chain amino acid degradation (bottom right). 



diguanylate cyclases (GGDEF domain) have been implicated as sig- 
nal molecules that trigger changes in the bacterial cell surface and 
motility 21, 27 . Some evidence suggests that cyclases also serve as sen- 
sors for osmotic pressure in eukaryotes 28 . Signal transduction 
domains, the weak homology of Sama2418 to a flagellar hook-assoc- 
iated protein, and the timing of expression of this regulator suggest 
that this protein may serve as an osmotic sensor in SB2B. In addition, 
protein Samal653 is a potential methyl-accepting chemotaxis sens- 
ory transducer protein with transmembrane characteristics, a signal 
peptide, a kinase domain, and a HAMP domain that ultimately 
results in the regulation of rotation of the flagellar motor 29 . The 
clustering of expression patterns of the regulator (Sama2418) and a 
methyltransferase (Sama0543) (Fig. 5) suggest osmotaxis may be 
driven by methylation of a methyl-acceptaing chemotaxis protein 30 . 

Because cell volume and ion regulation are not instantaneous 
processes, osmotic stress may damage cellular macromolecules and 
impair cell function until compensatory acclimations counteract the 
stress. Damage to DNA or changes in transcription caused by nega- 
tive supercoiling at high salt 31 33 leads to impairment of cell function 
and to the induction of repair processes and protection systems 34 . In 
SB2B, SbcC (Samal268), an exonuclease that removes unusual DNA 
structure(s), such as hairpins, that may be generated upon DNA 
damage 35 or through supercoiling is more abundant with salt stress. 
To reset the proteome to function in an environment of high salt, 
ClpA has been shown to have similar function to the DnaK and DnaJ 
chaperones during replication 36 . However, under osmotic stress 
induced by exposure to salt, the binding of ClpA to helicases involves 
a switch from a salt- sensitive to a salt-insensitive protein complex 37 . 
Although DnaK and DnaJ are often induced immediately following 
osmotic stress, our observations suggest that ClpA is associated with 
long-term acclimation and not the immediate stress response. 



As reported from a transcriptome analysis of S. oneidensis 9 , salt 
stress response and branched-chain amino acids (isoleucine, leucine, 
and valine) metabolism seem to be linked, prompting a hypothesized 
decrease in branched-chain fatty acid content at high salinity 38 as is 
found in other organisms 39 . We found that between 15 and 30 min- 
utes following salt exposure, SB2B expressed proteins that degraded 
branched chain amino acids completely to central metabolic path- 
ways (Fig. 3). However, we found no repression of any of the proteins 
involved in the synthesis of branched chain fatty acids that would be 
expected with a shift in the branched chain fatty acid content of the 
cell. Furthermore, fatty acid profiles of SB2B in response to high salt 
show no significant change in the concentration of branched chain 
fatty acids (Table 3), nor do original studies on the compositional 
changes of the fatty acid profile in Shewanellae during exposure to 
salt when grown on nutrient rich media 40 . Change in branched chain 
fatty acid profile during salt stress for S. putrefaciens only occurs with 
the removal of all exogenous sources of amino acids 41 suggesting a 
shift to branched chain amino acid metabolism in the media and at 
the expense of fatty acid content when necessary. Expression profile 
similarities illustrated in (Fig. 5) suggest that Sama3141, a putative 
adhesion protein, may be involved in the metabolic shift to degrade 
branched chain amino acids for Shewanella at high salt. 

Despite previous claims for Shewanella is , SB2B does not selec- 
tively degrade branched chain fatty acids as a response to high salt 
concentrations. The role of branched chain amino acid degradation 
in the adaptation of SB2B remains unclear, however the evidence 
presented here, in context with confirmed reports of changes in 
branched chain fatty acid content of Bacillus 39 raises the question 
whether membrane fatty acid changes in other organisms are a res- 
ponse to high salt or the result of a metabolic shift to branched chain 
amino acid degradation. The advantage of the time-course approach 
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Table 3 | Fatty acid composition of S. amazonensis SB2B during 
acclimation to osmotic stress. Values represent percent of total 
fatty acids. 

Time following NaCl pulse 





control 
0 min 


1 5 min 


30 min 


60 min 


90 min 


1 1 :0 iso* 


nd 


nd 


nd 


0.3 


0.3 


1 1:0 


nd 


nd 


nd 


0.3 


nd 


12:0 


3.5 


2.9 


3.1 


2.9 


2.6 


1 1:0 30H 


nd 


nd 


0.9 


0.9 


0.8 


13:0 iso* 


4.5 


3.8 


4.3 


4.8 


4.6 


13:0 


0.7 


0.7 


1.0 


1 .2 


1.1 


12:0 iso30H* 


nd 


nd 


nd 


0.4 


0.4 


12:0 3OH 


3.7 


3.0 


3.1 


2.8 


2.5 


14:0 iso* 


0.5 


0.5 


0.7 


1.1 


1.3 


14:0 


2.9 


2.8 


2.9 


2.5 


2.3 


13:0 iso30H* 


5.2 


4.3 


5.0 


5.5 


5.3 


13:0 301-1/15:1 isoH* 


0.4 


0.4 


0.6 


0.6 


0.6 


15:0 iso* 


25.7 


25.8 


25.5 


25.2 


25.9 


1 5:0 anteiso* 


0.8 


0.8 


0.8 


0.8 


0.8 


15:1 co8c 


nd 


0.3 


0.4 


0.4 


0.4 


15:1 co6c 


nd 


nd 


nd 


0.3 


0.3 


15:0 


2.1 


2.3 


2.4 


2.8 


3.2 


14:0 30H/1 6:1 iso 1* 


0.5 


0.4 


0.5 


0.4 


0.4 


16:0 iso* 


0.7 


0.8 


0.8 


1.2 


1.5 


16:1 cn7c/16:l G)6c 


19.9 


19.8 


17.1 


15.9 


15.5 


16:0 


12.1 


12.5 


11.1 


9.5 


9.4 


1 7:1 iso <b9c* 


nd 


0.4 


0.4 
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taken in this study is that it provides expression profiles that can be 
compared to determine co-expression patterns using multivariate 
statistical algorithms. The temporal expression patterns of the 
SB2B proteome following osmotic shock suggest an immediate arrest 
of transcription and division, followed by a shift to branched chain 
amino acid degradation, and the resumption of cell division. High 
expression of an adhesion protein that groups with branched chain 
amino acid degradation may suggest that the down-regulation of 
flagellar genes 9 and proteins reported here may coincide with a life- 
style change in Shewanella at high salt concentrations. 

Methods 

Media and growth conditions. Initial experiments were carried out on 
S. amazonensis SB2B (ATCC: BAA- 1098) to determine at what concentration of 
sodium chloride limits SB2B growth in LB broth. Initial cultures were inoculated to 
15 mL LB broth with NaCl concentrations ranging from 1% (0.17M NaCl; normal 
salt concentration of LB) to 8% (w/v; 1.37M) at 30 U C on an orbital shaker set at 
200 rpm. Growth rates were determined spectrophotometrically (OD 600 nm) in 
triplicate over a 24 hour period. We selected a NaCl concentration corresponding to a 
50% reduction in the growth rate for further proteomic analyses. 

To induce an osmotic stress response, cells were cultured in LB broth until log- 
arithmic phase growth {OD 600 nm —0.4) was achieved, followed by addition of NaCl 
to a final concentration of 0.85M salt (5%). The acclimation time was determined 
from the point of salt addition until resumption of growth. Growth rate dynamics 
were monitored by optical density and confirmed using standard plate count tech- 
niques. 

Sample preparation. Based on the acclimation time of SB2B to salt, samples were 
harvested at five time points following the addition of salt; 0 (control immediately 
prior to addition of NaCl), 15, 30, 60, and 90 min. Five replicate samples were 



harvested by centrifugation (5,000 rpm, 4°C, 10 min) and the cells were snap frozen 
using liquid nitrogen and stored at — 80 'C. 

For protein extraction, cell pellets were thawed on ice, washed, and suspended in 
100 raM NH4HCO3 (pH 8.4) then lysed via bead beating. From the cell lysis material, 
proteins were extracted and digested using trypsin into peptides according to estab- 
lished protocols 42 . Resulting peptide samples were pooled to create a composite 
sample for each time point to best utilize allotted instrument resources. 

Total RNA was extracted from cell pellets collected from the same time points 
using an RNeasy mini kit (Qiagen, Valencia, CA), and the RNA was quantified 
spectrophotometrically. Approximately 0.1 to 1 p.g of total RNA was incubated with 
genomic DNA elimination buffer and converted to cDNA using a QuantiTect reverse 
transcription kit (Qiagen, Valencia, CA). 

Proteomics data generation. Proteomics data was generated using the accurate mass 
and time (AMT) tag proteomics approach 16 . In brief, this approach utilizes tandem 
mass spectrometry to generate a reference peptide database (accurate mass and time 
tag database; AMT tag database) of observed peptides and their associated theoretical 
masses and elution times (normalized). This database is utilized to assign peptide 
sequences to ion current (relative abundance) information to peptides measured 
using high- resolution, high mass measurement accuracy mass spectrometry (LC- 
MS). The AMT tag database for S. amazonensis SB2B is maintained at the Pacific 
Northwest National Laboratory (Richland, WA) and was constructed as part of 
research efforts on Shewanella by the Shewanella Federation 2 ' 42 " 43 , For the study 
reported here, we utilized the normalized elution times and computed masses of 
peptides within the database to match sequences to peptides and their abundances 
measured using LC-MS from the NaCl perturbed and control culture samples. 

As part of this approach, samples for LC-MS analysis were first randomized and 
blocked; where, a single block contained one proteomics analysis per time point for a 
total of five blocks. Each block of samples was analyzed using an LTQ-Orbitrap 1M 
(ThermoFisher Scientific, San Jose, CA) mass spectrometer interfaced with a reverse 
phase HPLC system for peptide separation (LC-MS). Peptides were reverse-phase 
separated on in-house manufactured columns (60 cm X 360 um o.d. X 75 um i.d. 
fused silica capillary tubing) packed with 3 jim Jupiter C 18 stationary phase 
(Phenomenex, Torrence, CA). The HPLC system was equilibrated with 100% mobile 
phase A (0.2% acetic acid and 0.05% TFA in water). Fifty min. after peptide injection, 
mobile phase B (0.1% TFA in 90% acetonitrile/10% water) displaced mobile phase A, 
generating an approximate exponential gradient. Split flow controlled the gradient 
speed operating under constant pressure (10 kpsi). Separated peptides were ionized 
(positive) using an electrospray ionization (ESI) interface (manufactured in-house) 
with chemically etched electrospray emitters (150 mm o.d. X 20 mm i.d). 

Separated peptides were ionized (positive) using an electrospray ionization inter- 
face (manufactured in-house) that consisted of chemically etched electrospray 
emitters (150 mm o.d. 20 mm i.d). The LTQ-Orbitrap™ mass spectrometer was 
operated using a heated capillary temperature and spray voltage of 200°C and 2.2 kV, 
respectively. Data was acquired for 100 min, beginning 60 min after sample injection 
(10 min into gradient). LTQ-Orbitrap™ spectra (AGC 1X10 6 ) were collected from 
400-2000 m/z at a resolution of 100k followed by data dependant ion trap generation 
of MS/MS spectra (AGC 3X10 4 ) of the six most abundant ions using a collision 
energy of 35%. A dynamic exclusion time of 60 sec was used to discriminate against 
previously analyzed ions. Generated MS/MS spectra were searched using the 
SEQUEST algorithm 44 against the publicly available S. amazonensis SB2B translated 
genome sequence (www.jgi.doe.gov). Identified peptides of at least six amino acids in 
length, having 1 + , 2 + , and 3 + charge states and SEQUEST cross -correlation (Xcorr; 
an indicator of the statistical match between the observed and theoretically generated 
MS/MS spectra) scores of > 1.9, 2.2, or 3.5, respectively, were added to the existing 5. 
amazonensis SB2B AMT tag database. 

Orbitrap™ spectra were deisotoped using the software tool Decon2LS 45 and the 
tool VIPER 44 was used to identify mass and elution time features within mass mea- 
surement accuracy and elution time accuracy cut-offs of 5 ppm and 1%, respectively. 
Abundances for these features were calculated by integrating the signal strength 
under each peak of the LC-MS spectra 17 . To assign peptide sequences to this quant- 
itative data, mass and elution time features were matched to peptide sequences within 
the AMT tag database (identified by LC tandem mass spectrometry) for S. amazo- 
nensis SB2B. Peptide sequences were filtered according to a PeptideProphet™ 46 
probability value of < 0.90 prior to the matching of mass and elution time features. 
Matched features from each Orbitrap™ analysis (dataset) were then filtered on a false 
discovery rate (FDR) of less than or equal to 5%; the FDR associated with the AMT tag 
proteomics approach, calculated using STAC (Statistical Tools for AMT tag confid- 
ence). STAC is the latest statistical algorithm for assigning confidence to matched 
mass and elution time features, and is an evolution of an earlier algorithm referred to 
as SMART 47 . A protein was considered positively identified if a minimum of 2 unique 
peptides from its sequence were detected by LC-MS. 

Proteomics data analysis. The computational strategy set forth by Du, et al. 12 for 
analysis of label-free temporal proteomics data was followed. In brief, this strategy 
encompasses peptide filtering across replicates, abundance normalization, missing- 
value imputation, protein abundance estimation from peptides (protein rollup), and 
trend analysis. The associated MatLab™ (MathWorks, Inc., Natick, MA) code 
developed by Du, et al. (http://omics.pnl.gov/software/) was utilized 22 ; with the 
exception that central tendency normalization was performed using DANTE (http:// 
omics.pnl.gov/software/) 23 . A final trend for each protein was calculated as the 
median of trends observed from the replicate LC-MS measurements. 
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Gene expression analysis. Primers for genes involved in branched chain amino acid 
degradation, TCA cycle, motility, and DNA replication/repair were designed using 
the Primer 3 application (http://frodo.wi.mit.edu/primer3/) with a targeted amplicon 
size of 80-100 bp (Table 1). One ng ofcDNA was run in triplicate for 40 cycles, using a 
two step PCR in an Opticon 2 DNA engine (MJ Research, Waltham, MA). cDNA was 
labeled using a QuantiFast SYBR PCR kit (Qiagen, Valencia, CA) and 125 nM of each 
primer. Amplicon size and reaction specificity were confirmed by agarose gel 
electrophoresis and product dissociation curves. The number of target copies in each 
sample was interpolated from its detection threshold (C T ) value using a purified PCR 
product standard curve. The standard curve consisted of a serial dilution of a purified 
16S rRNA gene sequence from SB2B of known concentration. A representative 
amplicon for each gene tested was sequenced to confirm the target. 

Fatty acid analysis. In order to measure the membrane acclimation of SB2B to 
osmotic stress, cell samples used for gene expression and proteomic analyses were 
analyzed for fatty acid content. Triplicate samples collected for each time point were 
pooled and saponified using 45 g NaOH, 150 ml methanol and 150 ml water 
followed by methylation using 6.0 N hydrochloric acid and methanol. Fatty acids 
were extracted using a 1 : 1 solution of hexane and methyl tert-butyl ether and 
analyzed for composition as per the MIDI protocol (Microbial ID, Newark, DE). 

Motility. Cell motility was qualitatively tested in triplicate using soft agar 
inoculations on LB plates. We prepared solid (1% agar) and soft (0.3% agar) plates 
both at salt concentrations of 0.17M (normal LB media salt concentration) and 0.5M 
for motility assessments. Each plate was inoculated with SB2B (5ul; OD 600 — 0.5) and 
incubated at 25°C for 72 h. Colony diameter was measured for comparison. 
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